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" ^ , The title 'Chiral extrapolations' of this contribution refers more precisely to chiral ex- 
trapolations of lattice QCD data. We shall deal with low-energy aspects of QCD, which 
^—1 . are not accessible to ordinary weak-coupling perturbation theory. Possible alternatives to 
■ weak-coupling perturbation theory in the low-energy domain of QCD include the investi- 
gation of specific models, Monte Carlo simulations of lattice regularised QCD, and chiral 
effective field theories (CHEFT), the latter being low-energy theories which incorporate 
(«~») the constraints from (spontaneously broken) chiral symmetry. It is the comparison of 
"j—^ CHEFT (or chiral perturbation theory (CHPT)) with results from lattice QCD simulations 
d that will be the subject of the present paper. 

' However, the reader should be warned that as lattice QCD practitioners we look at 

CHEFT from an (ab)user's viewpoint. Also a second warning may be in order: This is 
(~| not a review. The material has been selected according to subjective criteria, and the 
references are far from being complete. For a review of CHPT see, e.g. , (Leutwyler 2000) 
. ^ and (MeiBner 2000). For reviews on closely related subjects, partly overlapping with the 
present work, see (Bar 2004) and (Colangelo 2004). 

_ _ I 

2 Lattice regularisation and Monte Carlo simulation 



The basic input for a Monte Carlo simulation of lattice QCD is first of all a lattice action, 
i.e. a discretised version of Euclidean QCD. Secondly, one has to choose a lattice size 
(necessarily finite). Thirdly, the bare coupling constant and the quark mass(es) have to 
be fixed. Then the lattice spacing a and the spatial box size L can (approximately) be 
given in physical units. 

In general it is not possible to choose the simulation parameters such that the results 
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can immediately be identified with experimentally measurable quantities. In particular, 
three extrapolations are required: the continuum limit a ^ 0, the thermodynamic limit 
L 00, and the chiral hmit, where the masses of the light quarks decrease to their physi- 
cal values and further down to zero. Unfortunately, in all three limits the simulation costs 
increase rapidly. It is therefore preferable to appeal to theory in order to relate simula- 
tion results obtained for unphysical quark masses, finite volumes, ... to phenomenology. 
This wiU then lead to well-justified extrapolation formulae. 



3 Chiral effective field theory in the pion sector 

The natural starting point for chiral perturbation theory is the pion sector. The very exis- 
tence of Ught pions (for the case of two flavours) relies on the spontaneous breakdown of 
chiral symmetry combined with the weak explicit breaking due to the quark masses. A 
further consequence is the weakness of the pion-pion interaction at low energies and mo- 
menta which makes a perturbative treatment meaningful. This (chiral) perturbation the- 
ory is most conveniently set up by means of an effective Lagrangian, i.e. the most general 
Lagrangian for effective pion (and later on also nucleon, . . . ) fields which is compatible 
with chiral symmetry. The effective Lagrangian is constructed out of terms with more 
and more derivatives as the order of the expansion increases, leading at the same time to 
an increasing number of effective coupling constants, which are not determined by chiral 
symmetry. Eventually one obtains an expansion of the physical observables in the pion 
mass and the particle momenta, where both are considered as quantities of the order 
of a small parameter p. Here 'small' means p <C 47rF^ ^ 1 GeV. 

We shall restrict ourselves to the case of two flavours, Nj = 2, with isospin breaking 
neglected, i.e. we take for the quark masses m„ = = niq. Then the lowest-order ex- 
pression for is the famous Gell-Mann-Oakes-Rennerrelation = 2\{qq)\mq/F'^, 
where (qq) is the chiral condensate and F denotes the pion decay constant in the chi- 
ral limit. 

Beyond leading order one finds (see, e.g., Colangelo et al. 2001) 

M2 = m2 {l - \xh + fx^ili + x^kM + 0(x')} (1) 

with 

Im = ^ {28h + 324 - 94 + 49) . (2) 
51 

Here we have set M2 = 2mgBwithB = - {qq)/F'^ and a; = M'^ / [IQ-k"^ F"^) . The chiral 
logarithms are hidden in the quantities £i = ln(Af /M^), which contain the information 
on the (renormalised) coupling constants. The term proportional to Rm represents ana- 
lytic contributions 0{x'^), which are expected to be small. Note that Um and Aj depend 
neither on niq nor on the renormalisation scale. One can estimate F = 0.0862 GeV, 
and phenomenological analyses lead to Ai = 0.121°;°^ GeV, A2 = 1.20t[J [J^ GeV, 
A3 = 0.591 J:^? GeV, A4 = 1.25lg:i^ GeV (see, e.g., Colangelo and Durr (2004)). 
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4 Comparison with Monte Carlo data in the pion sector 



Let us start with some general remarks on the comparison of Monte Carlo data with 
CHEFT formulae. First of all, one needs results in physical units. A popular way to set 
the physical scale uses the Sommer parameter rg, which is a length scale derived from 
the heavy-quark potential V{R) through the condition dV{R)/dR\ii=ro = 1-65. The 
phenomenological value has been found to be approximately tq = 0.5 fm, which is the 
number to be used in the following. Note that this method assumes that the dependence 
of ro on the light quark masses is negligible, an assumption whose validity is not quite 
clear Secondly, many quantities, such as, e.g., quark masses, have to be renormalised. 
Finally, we have to deal with lattice artefacts. Ideally, one would eliminate them by an 
extrapolation to the continuum limit, which is not an easy task, or one could incorporate 
them in the CHEFT (see the review by Bar (2004)). In the following we shall adopt a 
simple-minded approach and try to select the data such that cut-off effects are negligible. 

Typically there is little structure in the quark-mass dependence of lattice results, and 
the data are in most cases well described by a linear function of niq. In other words, 
there are no obvious chiral logarithms. Of course, this may be due to the relatively large 
quark masses in the present simulations, where leading order CHPT is unlikely to work. 
Thus one needs higher-order calculations, and one has to face the question up to which 
masses CHPT is reliable. Alternatively, one may try to tame the unphysical behaviour of 
the truncated series at large masses by some cut-off function and in this way arrive at a 
formula which works in the mass range covered by the simulations. For this approach 
see, e.g., (Leinweber et al. 2004), (Young et al. 2004), and references therein. 

In Fig. [Owe compare pion and quark masses obtained by the JLQCD collaboration 
(Aoki et al. 2003) with the quark-mass dependence of the pion mass as predicted by 
Eq. The parameters have been chosen (not fitted) as follows: F = 0.0862 GeV, 
B = 3.8 GeV, Ai = 0.12 GeV, A2 = 1.20 GeV, A3 = 0.65 GeV, A4 = 1.25 GeV. It is 
a remarkable observation that the Gell-Mann-Oakes-Renner relation — 2Bniq is a 
rather good approximation for pion masses up to about 0.7 GeV. The above parameters 
are well compatible with the phenomenological values quoted in the preceding section. 
Only B is somewhat larger than the value B — 2.8 GeV given in (Diirr 2003). Note 
that B is scale and scheme dependent as are the quark masses, while the product Bniq 
is independent of scale and scheme. Here we have employed tadpole improved one- 
loop perturbation theory with the renormalisation scale fi — 2 GeV in order to convert 
the bare VWI masses to renormalised quark masses niq in the MS scheme. The use of 
perturbation theory entails a considerable uncertainty in the renormalised quark masses 
and hence in B. Indeed, a recent investigation (Gockeler et al. 2004) suggests that the 
non-perturbative mass renormalisation factor (at the bare coupling used by the JLQCD 
collaboration) is about 2.3 times larger than the perturbative estimate employed here. But 
it is gratifying to see that chiral perturbation theory with phenomenologically acceptable 
values of the coupling constants is able to make contact with the low-mass end of the 
quark mass range that can be reached in present simulations with dynamical quarks. For 
a more detailed discussion of the quark-mass dependence of in comparison with 
different Monte Carlo data see (Diirr 2003). 

For the pion decay constant (normalised such that the physical value is 92.4 MeV) 
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Figure 1. The square of the pion mass versus the quark mass. JLQCD data for Nf = 2 
are compared with the Gell-Mann-Oakes-Renner relation (dotted line) and chiral per- 
turbation theory (full line). 



chiral perturbation theory yields (see, e.g., Colangelo et al. 2001) 

F^ = Fi^\ + xh - \x^ij, + x^kp + 0{x^)^ (3) 

with 

ip = ^{lAh + IQh + 64 - 6I4 + 23) = In ^ . (4) 

Again, the analytic contributions 0{x^) (proportional to kp) are expected to be small. 
In Fig. 12] we compare this formula with preliminary data from the UKQCD and QCDSF 
collaborations. Motivated by our observation that the Gell-Mann-Oakes-Renner relation 
works so well, we replace M by M^. Thus we avoid the problem of renormalising the 
quark mass. Choosing F = 0.0862 GeV, kp = -0.5, A4 = 1.12 GeV, Af = 0.67 GeV 
(consistent with phenomenology) we obtain the full curve in Fig. |2 which connects the 
physical point with the data for the lowest mass but does not describe the data at larger 
masses. It is however reassuring that for masses up to the first Monte Carlo points the 
chiral expansion seems to be well-convergent: The dotted curve, which corresponds to 
~ F{1 + xii}, does not deviate dramatically from the full curve in this region. 
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Figure 2. Preliminary data ( Nf — 2) for the pion decay constant compared with chiral 
perturbation theory at different orders. The asterisk indicates the physical point. 

5 Including baryons 



The nucleon mass mjv does not vanish even in the chiral Hmit. Indeed 771 tv ^ M^, and 
a non-relativistic treatment of the nucleon field seems reasonable. This leads to the so- 
called heavy-baryon chiral perturbation theory (hbchpt). A relativistic formulation of 
chiral perturbation theory in the baryon sector has been given by Becher and Leutwyler 
(1999). In both cases the physical picture is that of a (heavy) nucleon core surrounded by 
a cloud of light pions. This is rather different from the situation for the pion, and hence 
the behaviour of the chiral expansion in the nucleon sector need not be similar to that in 
the pion sector. 

In Becher's and Leutwyler's formulation one obtains for the nucleon mass (Becher 
and Leutwyler 1999, Procura et al. 2004) 



mjy = niQ — AciM^ 
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(5) 

Here we have again identified Af2 = M'^ = 2Bmq, A is the renormalisation scale, mo 
and are the mass and the axial charge of the nucleon in the chiral limit, ci, C2, C3 
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Figure 3. Nucleon mass data for Nf — 2 on (relatively) large and fine lattices. The 
asterisk indicates the physical point. The curve corresponds to the fit mentioned in the 
text. 



denote coupling constants from the effective Lagrangian, and e\{\) is a counterterm. 

Hadron masses for Nf — 2 have been published by the CP-PACS collaboration 
(Ali Khan et al. 2002), the JLQCD collaboration (Aoki et al. 2003) as well as the UKQCD 
and QCDSF collaborations (Allton et al. 2002, Ali Khan et al. 2004a). We have selected 
results obtained on (relatively) large and fine lattices to compare with chiral perturbation 
theory. More precisely, we have considered masses from simulations with a < 0.15 fm, 
< 800 MeV and Mj^L > 5. These ten data points were fitted with Eq. (|3 where 
g^, F, C2 and C3 were fixed to phenomenologically reasonable values while mo, ci and 
e\{X = IGeV) were the fit parameters. For more details see (Ali Khan et al. 2004a). 
The fit curve and the data points are shown in Fig.|3l It is first of all remarkable that the 
masses obtained by different collaborations with different lattice actions and algorithms 
fall (to rather good accuracy) onto a single curve. Furthermore, the fit parameters are 
very well compatible with phenomenology, in particular, ci is about —1 GeV^^ and the 
fit curve comes quite close to the physical point. On the other hand, Eq. ^ seems to 
work up to surprisingly large masses. 



6 Axial charge of the nucleon 

The quark-mass dependence of the axial charge (or axial-vector coupling constant) 
of the nucleon has been studied by Hemmert et al. (2003) within the framework of the 
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so-called small-scale expansion. This is an extension of HBCHPT which includes explicit 
A(1232) degrees of freedom. 

In the small-scale expansion the expansion parameter is usually called e, and one 
finds in O(e^): 



(g^A^M^ 
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The new parameters appearing here are (the value of gA in the chiral limit), A (the 
nucleon A mass splitting in the chiral limit), ca, 5i (the A^A and AA axial coupling 
constants), and C^^^{X) (a counterterm at the renormalisation scale A). Hemmert et al. 
(2003) find that for reasonable values of the parameters the formula (|6j is able to describe 
the rather weak mass dependence of the Monte Carlo data as well as the physical point. 



7 Chiral effective field theory in a finite volume 

Presently, lattice QCD simulations are not only restricted to unphysical quark masses, 
but also to relatively small (spatial) volumes, usually with periodic boundary conditions. 
While simulations at the physical quark masses might be possible some day, it will take 
a bit longer before the ideal case of an infinite volume simulation can be realised. In 
the meantime we can take advantage of the fact that the chiral effective Lagrangian is 
volume independent for periodic boundary conditions (Gasser and Leutwyler 1988). So 
the same Lagrangian governs the quark-mass as well as the volume dependence, and 
additional information on the coupling constants can be extracted from finite size effects. 
This description of the finite size effects should work as long as they result from the 
deformation of the pion cloud in the finite volume, i.e. as long as L is not too small. 
After all, it is the pion propagation that is predominantly affected by the finite volume, 
because the pion is the lightest particle in the theory. Treating L^^ as a quantity of order 
p Uke Mjr we arrive at the so-called p expansion (Gasser and Leutwyler 1987). This is 
to be distinguished from the e expansion, where L^^ — 0(e), M-^ = O(e^) with a small 
parameter e. 

In more technical terms, the finite volume (with periodic boundary conditions) discre- 
tises the allowed momenta such that the momentum components are restricted to integer 
multiples of 2tt/L. So the loop integrals of CHPT become sums. On the other hand, we 
can interpret the resulting expressions in the following way. In a finite volume a pion 
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Figure 4. Volume dependence of the nucleon mass for M-^ = 545 MeV. The dotted 
curve shows the contribution of the term, while the solid curve includes also the p'^ 
correction. 



emitted from a nucleon can not only be reabsorbed by the same nucleon, but also by one 
of the periodic images of the original nucleon at a distance which is an integer multiple 
of L in each of the finite directions. From this point of view the finite size effects arise 
from pions travelling around the volume once, twice, ... in a given direction, and each 
crossing of the boundary leads (roughly) to a factor exp{—M.^L) in the final contribution. 

As an example, let us consider the nucleon mass in relativistic baryon CHPT. Then 
the nucleon mass in a spatial box of length L, but for an infinite extent in time direction, 
can be written as 



where 



mjv(i) = mAr(oo) + A„(L) + A,(L) + 0(p5), (8) 



AJL) ^^-^M^^Mlj^A^ Y: [L\n\^ml.^+Ml{l-x)^ (9) 



is the 0{p^) contribution to the finite size effect and 



Ki{L\n\A'U) , K2{L\n\M^) 



(2C1 - C3) rr^TTT h C2 



L\n\M^ {L\n\M^y 



(10) 



is the additional contribution arising at 0{p'^). Here can be interpreted as the number 
of times the pion crosses the 'boundary' of the box in the i direction. Note that the 
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finite volume does not introduce any new coupling constant. For a comparison with 
Monte Carlo data we evaluate the finite size corrections using the parameters from the 
fit discussed in Section|5] choose itin {oo) such that mN{L) on the largest lattice agrees 
with the Monte Carlo value and take for A/tt the value from the largest lattice. This 
leads to a surprisingly good agreement with the Monte Carlo masses, as is shown for 
= 545 MeV in Fig.0]using JLQCD data. Even for larger pion masses the formula 
reproduces the finite size effects quite well (Ali Khan et al. 2004a). This lends further 
support to the fit shown in Fig.^] 

In the case of the pion mass the leading finite volume correction has already been 
known for some time (Gasser and Leutwyler 1987): 



Moreover, neglecting pions which travel around the box more than once (|n| > 1) 
Liischer (1986) has expressed the finite volume correction in terms of the tttt forward 
scattering amplitude F{i'), i/ = (s — u)/{4Mt^): 



The chiral expansion of F{i/) is known to 0{M^) and the corresponding coupling con- 
stants are reasonably well determined. Using this information one can work out the 
volume dependence of Af^ predicted by Eq. (I12> (Colangelo and Diirr 2004). The non- 
leading terms of the chiral expansion give non-negligible contributions, still the finite size 
effects are considerably underestimated by this approach, at least for A/jr > 500 MeV. 
Probably terms with \n\ > 1 are important for smaller volumes, which is certainly the 
case in Eq. (^3- However, Eq. (II 1> alone predicts even smaller finite size effects than 
Eq. (I12t . So it seems that one would need higher orders in the chiral expansion as well 
as pions propagating around the volume more than once, and the present understanding 
of the finite size effects for A/^r is unsatisfactory. For another investigation of hadron 
masses in a finite volume see (Orth et al. 2004). 

The volume dependence of has recently attracted much interest. First calculations 
in chiral perturbation theory have been performed by Beane and Savage (2004). How- 
ever, at the masses used in current simulations these leading-order formulae do not even 
reproduce the sign of the finite size effects observed in the Monte Carlo data, see, e.g., 
(Ali Khan et al. 2004b), (Sasaki et al. 2003). It remains to be seen whether more ad- 
vanced calculations in CHEFT can solve this discrepancy or whether lower quark masses 
are required. 

8 Conclusions 

What has been described here, are only the first steps of an ongoing effort to combine 
calculations in CHEFT and lattice simulations. The overall impression is that the range 
of quark masses that can be used in actual Monte Carlo computations is beginning to 
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overlap with the region of applicability of CHPT. In favourable cases this overlap seems 
to be so large that meaningful fits are possible. Although fits without phenomenological 
input are still beyond reach, interesting information on (effective) coupling constants can 
already be extracted. In most cases the leading chiral logarithm appears to be dominating 
only for rather small quark masses. Thus CHEFT should be pushed to higher orders while 
on the lattice side results for lower quark masses are eagerly awaited. 
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